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Cities are living organisms. They are out of equilibrium, open systems that never stop developing 
and sometimes die. The local geography can be compared to a shell constraining its development. In 
brief, a city's current layout is a step in a running morphogenesis process. Thus cities display a huge 
diversity of shapes and none of traditional models from random graphs, complex networks theory 
or stochastic geometry takes into account geometrical, functional and dynamical aspects of a city in 
the same framework. We present here a global mathematical model dedicated to cities that permits 
describing, manipulating and explaining cities' overall shape and layout of their street systems. 
This street-based framework conciliates the topological and geometrical sides of the problem. From 
the static analysis of several French towns (topology of first and second order, anisotropy, streets 
scaling) we make the hypothesis that the development of a city follows a logic of division / extension 
of space. We propose a dynamical model that mimics this logic and which from simple general rules 
and a few parameters succeeds in generating a large diversity of cities and in reproducing the general 
features the static analysis has pointed out. 



I. INTRODUCTION 

The city is a living structure: it is an open system, 
always in motion. A city is born, develops, heals over 
injuries (war damages...) and sometimes dies in part or 
totally. Its development responds to internal and exter- 
nal constraints in such a way that local geography acts as 
a shell that sculpts it. As general living systems, cities 
exhibit a huge range of diversity both on their overall 
shape (that can be circular, sprawling, linear or even 
fractal) and on the appearance of their street systems 
(regular, organic, tree-like). Such a diversity can also be 
observed at the level of a single city that has not devel- 
oped homogeneously. 

We seek out to show that behind this diversity stands a 
single principle: a city develops within a logic of division 
/ extension of space. The intra and inner diversity of 
shapes can be seen as a variation in a coherent global 
phenomenon. Our approach is street-based: we consider 
an infinitesimal piece of street as the elementary compo- 
nent of a city and bet it contains important information. 
The dual approach is to consider the build-up area (build- 
ings, parks...) as the unit of formation. Several points 
of view have been used in the past to model cities: cel- 
lular automata, multi-agent systems, fractals, stochastic 
geometry, L-systems and graph theory leading to com- 
plex systems' theory. 

Cellular automata and multi-agents systems have widely 
and successfully been used to simulate the dynamics of 
populations and of land use [Tj. The fractal description 
of cities [2] gathers these simulations into a theory and 
point out the advantages of a fractal shaped city from 
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the point of view of the build-up area. Nonetheless the 
basis of these models is either a discrete field or the map 
of a given city. They explain the global differentiation of 
space when the street network is known or ignored. 
Since the famous "Bridges of Konigsberg" problem by 
Euler, one is tempted to describe a city as a graph with 
streets as edges and their intersections as vertices. This 
provides a relational representation of the city [3] . One of 
the difficulties then is the particular embedding of these 
graphs that make random graphs unable to stick to a 
city's map representation. 

Stochastic geometry gets around this problem by consid- 
ering stationary tessellations (Poisson Voronoi, Poisson 
Delaunay or Poisson Line Tessellation, their superposi- 
tions and iterations) that are geometrical objects em- 
bedded in a compact subset of the space and deduces 
geometrical random graphs from their induced topology 
[3] • L-systems with procedural programming make a map 
evolve from local coherence rules and input data that in- 
corporate global constraints 5,6. The stochastic geome- 
try approach gets good results at analyzing optimization 
problems on street networks and L-system are success- 
fully used in graphics but they do not explain the under- 
lying phenomena at work to determine the appearance of 
a city. 

For a few years a complex networks based study has been 
adopted to describe cities [7J [5] . But the main conclu- 
sion is that a city behaves neither like a classical scale-free 
network nor a small-world essentially because of its spa- 
tial embedding [5]. Cities then clearly need a dedicated 
mathematical framework taking into account both the 
topology and, what is new, the geometry of a street net- 
work. The scope of this article is a street-based approach 
of cities that allows analyzing, manipulating and explain- 
ing cities' morphogenesis. In pl| we define a mathemat- 
ical formalism to handle with cities both on a relational 
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and a geometrical way. The very stuff of the geometry 
being that street segments are lined up into sets called 
stree ts we model via a hypergraph additional structure. 
( |III[ ) presents several description measurements (topol- 
ogy of first and second order, anisotropy and streets scal- 
ing) to obtain quantitative comparison elements between 
cities and exhibit some features of the global mechanic 
of cities system (organic shape, log scaling of length and 
small topological radius). These features call for the 
modelling of the city as a process of division and exten- 
sion of space. We use 10 French towns and their centers 
without having to restrict them to a square window to 
illustrate this. 

Eventually we present a morphogenetic model of the city 
( IV I and its simulation ([v]) that implements the idea of 
space division / extension. This model reproduces the 



general features pointed out in (III) and the variation 



of a few normalized parameters allows recovering a large 
range of diversity. The model is expressed in terms as 
relevant for mathematicians or physicists as for town- 
planners or social scientists. No sampling of space is 
required for the simulation. 



II. CITIES' SPACE 

Since the famous resolution of the " Bridges of Konigs- 
berg" problem by Euler [3 one is tempted to look at 
cities with a formal and relational point of view: a city 
is a graph whose edges are streets and vertices are their 
intersections. 

Nevertheless this approach does not take into account 
the physical constraints (of being a functional object on 
the plan) that are exerted on cities. Furthermore, it en- 
hances a fundamental disjunction between intersections 
- that would be objects of interest - and street segments 
that would simply bind them. Under those street seg- 
ments / edges lies a characteristic geometry. Each point 
of this geometry should be seen as an object in relation 
with other similar objects. This paragraph aims at intro- 
ducing a terminological frame that permits manipulating 
cities as "continuous graphs" embedded in a two dimen- 
sional Euclidian space. The vocabulary used in this ar- 
ticle is freely adapted from general graph theory |10j to 
respond our to specific needs. 

When importing a map (via a .MIF file), the raw data is 
coded into a list of polylines. In fact these polylines come 
from the sampling of curved streets that are difficult to 
represent in a computer. We consider it is possible to 
transform a geometrical graph into a straight graph ar- 
bitrary close to it. The important paradigm of degree 
two points comes up. To rectify the map we have added 
degree 2 vertices but we want to consider a version the 
same object. In short, a segment with its two extremi- 
ties and the same segment with its extremities plus its 
midpoint should be seen as the same entity. We define a 
measure on a geometrical graph which allows to see that 
graph as a "continuum" and to consider each point of its 



geometry in the same time. The data rather represents 
street segments that actual streets. The particularity of 
a city's geometry is that its street segments are coher- 
ently arranged into disjoint geometrical sets: the streets. 
We try to record the notion of streets with a hypergraph 
additional structure. This provides a multi-scale repre- 
sentation of the city. 



A. Graphs and planar graphs 

Let S be a set, V (for vertices) a finite subset of S, E 
(for edges) a symmetric part of V x V then G = (V, E) 
is said to be a (undirected) graph. 

A drawing of G is an injective function from V to R 2 
and from E to the set of continuous paths such that the 
image of an edge has for limits the images of the vertices 
it binds and does not pass through images of vertices it 
does not bind. An edge-crossing is the intersection of the 
images of two edges outside the image of V. 
If there exists a drawing without any edge-crossing, the 
graph is said to be planar (Fig. [T|. The first character- 
istic of city graphs is their planarity. 



B. Geometrical graphs 

A geometrical graph can be seen as a particular draw- 
ing of a planar graph. 

Let the available space A be a connected and compact 
subset of R 2 , V a, finite subset of A and E a set of almost 
everywhere derivable paths included in A from one ele- 
ment of V to another that do not intersect outside of V. 
Then G = (V, E) is an element of the space of geometri- 
cal graphs <? 9 (A). If E is restricted to straight segments 
(G € (? S (A)) , G is a straight graph. To a geometrical 





FIG. 1. The representation of a graph by its symmetric 
adjacency matrix (a). This graph is planar: it admits at least 
two geometrical graphs as drawings (b). 
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graph G, one associates ttq, the subset of A defined by : 

7r G = {x e A, 3e e E, x e e} (1) 

-kg is compact so we can provide Q g (A) with an Hausdorff 
distance : 



<Msi || 52) 



max mm 



(2) 



A drawing G' = (V',E') is a rectification of the geo- 
metrical graph G = (V, E) if V C V C ttg and if each 
element of -B' is a segment. G' is not necessarily a planar 
straight graph since edges can possibly intersect outside 
of vertices (Fig. [2]). The idea is that one should be able 
to add to a geometrical graph as many vertices of degree 
two as he wishes and still consider the same mathemati- 
cal object. In this article we will admit that Prl and Pr2 
hold for a large enough class of geometrical graphs: 

Prl: : Every geometrical graph admits a planar rectifi- 
cation 

Pr2: : Every geometrical graph is the limit of a sequence 
of straight graphs 

To get (Prl), one has to sample an original graph with 
enough additional vertices of degree 2 and a small enough 
edges length. (Pr2) states that it is possible to approx- 
imate any geometrical graph by a straight graph arbi- 
trary close to it in the sense of dn • This is of practical 
importance as it allows to work only with straight planar 
graphs, while original city's maps can have curved edges. 





FIG. 2. A geometrical graph (a) and two of its rectifications 
(b). The upper one is not straight because two edges intersect 
outside of the vertices set. 



sure hq. For instance J G due is the total length of edges 
in G or in ( III C ) we define the angular density of a graph 
^(.) with Hg{ )- This measure respects: 



f f(g)d» G = J2\ I /W^i^l 



(3) 



Where fii is the i-dimensional borelian measure, / a pos- 
itive function and J the integral along the path e. If / 
is a continuous function defined on R 2 we have : 

/ fid) dfi G = lim ^ J f{x,y).{-K G ®B e )dn 2 {x,y) (4) 

where © is the Minkowski addition and B e the closed 
Euclidian ball of radius e: x © B e — { y g R 2 ||a; — 
y\\ < e }• So fie is a measure between fii and ^2 that 
allows to make quantitative measurements on the whole 
geometrical graph. 



D. Hypergraph structure 

If H is an equivalence relationship on E then 
((V, E),H) is said to be an hypergraph. 
Let (V, E) a graph and R a reflexive relationship on E 2 . 
Then the relationship R defined by : 

ei i?e 2 hf 3 a± = ei, a 2 , ■■■ , a n = e% G E | 

ati R a 2 , a 2 R a 3 , a n -i R a n (5) 

is an equivalence relationship. From this, one can con- 
sider Rg : 

ei R g e 2 iif (ei* 2 e 2 ) V((ei*e 2 ) A (|^(ei, e 2 )-7r| < 6)) 

where e\ *e 2 means that e\ and e 2 intersect, and e\ * 2 e 2 
that ei and e 2 intersect in a vertex of degree 2. Z(ei, e 2 ) 
stands for the angle between ei and e 2 oriented with the 
same origin. 

For instance think of considering the map of a city and 
the relationship "these two edges are pieces of the same 
street" . This Rg allows recovering the notion of "streets" 
even if input data do not contain such labels. The algo- 
rithm labeling streets segments with a street number does 
not depend on its starting point and is fast to run. The 
price to pay is that some special cases as forks of two 
segments making a very small angle with a third one will 
be considered as a single street. This over structure is 
essential as it gives a way to analyze overall structures of 
planar graphs, and in particular of cities. 



E. Citiy graphs 



Measure 



As a compact part of A, n G is a polish space (complete 
and separable) on which one can define a borelian mea- 



We define pragmatically the set of city graphs Q G as 
the subset of Qs that represents an existing city or a city 
that could have existed. We restrict this definition to Q s 
since a lot of cities' streets are straight and even if it is 
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FIG. 3. A straight graph (a) and its hypergraph structure 
(b) deduced from R^no- Viewed as a city's map, this graph 
contains 7 streets segments but 3 streets. 



not the case there exists a straight approximation as ac- 
curate as we want. 

In the following, we will write C — (V, E) a city and 
canonically provide C with an hypergraph structure from 
a relationship Rg. For the sake of simplicity we keep 
the same notation to designate the hypergraph : C — 
((V,E),H). Its borelian measure is written ^c- 
The relational aspect inherited from a simple graph struc- 
ture allows us to define the set of faces F. Euler's equality 
is respected so that: ((V — §E + $F = 1. To each edge 
e we associate the set V(e) of its extremities in V and 
to each point v £ nc we associate E{v) the set of edges 
that pass through it. N{v) is the degree of a vertex. 
We partition V = V\ U V% U V + where V\ contains all ver- 
tices of degree 1, Vi of degree 2 and V+ of higher degrees. 
Vertices in Vx, Vz and V+ are respectively terminations, 
junctions and intersections. Elements of V% will be seen 
has sampling artifacts used to fit the curved geometry of 
the city. 

Elements of E are called streets segments, those of H 
streets. 



III. 



CITY SPACE FEATURES 



Due to their spatial constraints, real geometrical 
graphs do not behave as classical complex networks 
(small- worlds or scale- free networks) j9]. Real cities dis- 
play structural, geometrical and functional features. For 
instance a real city aims both at lodging its inhabitants 
and at providing them with an efficient access to geo- 
graphical and human resources. These constraints log- 
ically affect the structure of a city graph. The pur- 
pose of this part is to define some mathematical tools 
that will quantitatively measure structural differences 
between city graphs. Classical measures from complex 
networks theory (efficiency, robustness, centrality, degree 
correlation) have been investigated in [TlTfl"3"] . The mea- 



sures presented below are dedicated to cities. Each one is 
illustrated with the city of Amiens in France: the whole 
city and its center (Fig. [5]) . The main properties are reca- 
pitulated in (HIE) for 10 French towns and their centers. 



We have imported vector maps in a calculus framework 
and successively: rectified them by taking care of con- 
serving the planarity and angles at the intersections and 
used -Rtt/io to obtain an hypergraph structure. 

A. First order topology 

Let C = ((V,E),H)) be a city and N(k) = %{v € 
C, N(v) = k} be the number of vertices of degree k in 
C and N(k) = N{k)/J2N(i). The set V 2 (junctions) 
should not be taken into account since it only represents 
sampling artifacts to preserve the shape of streets. In [8] 
the histogram of N is studied by means of an exponential 
tail of distribution. Nonetheless, this distribution (Fig. 
|4| is very peaked in 3 or in 4. It is sufficient to describe 
the histogram N by the organic ratio: 



JV(l) + iV(3) 



(7) 



which allows to discriminate quickly whether the city had 
been planned (r/v ~ in the limit case) or not. Indeed 
a planned city is filled with a homotopy of a rectangular 
grid (only N(i) ^ 0). This is clearly useful to settle 
buildings but also sticks to human perception of space 
since we have the intuition of left - right / front - behind. 
In unplanned cities (organic will be the dedicated word 
in IV to emphasize the comparison with living systems) 
i.e. created by the interaction between non concerted 
settlements, there is little probability for streets segments 
to be coherent and thus r^r — 1. 

Following [13] we characterize the topology of a city by 
its " meshedness coefficient". It's easy to count v and e 
then / is deduced from Euler's formula. Given V, the 
maximum number 2. s — 5 of faces is obtained by the 
Greedy Triangulation algorithm [13] . So the quantity 
M = (e — v + l)/(2.« — 5) equals to if the city is a tree 
and is close to 1 if it is a highly connected graph. For 
real cities [12] M typically ranges between 0.08 and 0.35. 
To be coherent with our preceding remark we should not 
take junctions into account. It has no incidence on the 
numerator but we have to change 2.v — 5 into 2.v.(l — 
N(2)) - 5 



1 



2.v.(l - N(2)) 



(8) 



M3 is quite small because of the general lack of triangles 
in the topology of a city. As a trapezoid contains two tri- 
angles, we rescale it in M 4 = 2M 3 whose maximum is hit 
when the considered city contains the maximal number 
of trapezoids. Amiens appears as an "average" organic 
city, with a meshedness coefficient of 0.41 (0.54 when re- 
stricted to the center) and — 0.79 (0.68 in the center). 
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FIG. 4. The histograms of degrees' distribution in Amiens 
(whole city and center). We observe that both distributions 
are peaked in 3 even in the more regular city centre, where 
the number of connections 4 is still half. The large number of 
2 comes from straightening, especially in the suburb curved 
streets and should not be taken into account. 



B. Second-order topology 

The streets H induce a particular topology. In [12] its 
study is referred as the " dual approach" . We use here the 
expression "second order topology" since it is a derived 
topology. Moreover, as faces are the dual of vertices, 
the dual of a space is a space containing the same infor- 
mation. We prefer to keep back the word "duality" to 
express this idea: " the mass of a city (buildings, houses, 
parks) is the dual of the street system" which could refer 
to the work of [DH]. 

Let us call topological distance of C the function d^ po : 
H x H — > N that satisfies: 



4 opo (/ii,7i 2 ) = 
d t t S, po {hi 1 h 2 ) = 



min dJ?°(h,h\) 

h£H,hnh 2 ^ftl 



if/ii = li2 
1 otherwise 

The topological distance counts the number of times one 
needs to turn to go from a street to another one. The 
topological average distance of a street ho is : 



(10) 



h£H 



This formula defines a new centrality measure on the map 
similarly to those studied in [11] . A street that minimizes 
dc P ° is then called a center. One drawback of this wholly 
topological definition is that, since topological distances 
are integers, several streets can be defined simultaneously 
as central streets. Common sense would then be to take 
all of these streets as simultaneous central streets, and 
to calculate the distance of any street as the minimum of 
the distance to any of these streets. Another way is to 
weight by the streets length: 



^topo leng 



(ha 



1 



\\h\\.d topo (h,h 



(ii) 



is the length of the street h). From this a unique 
center h c is defined if the city is not too regular. The 
topological radius of the city is then defined by r£ po = 
maxd^ pc '(h,h c ) and its diameter by : 



dtopo 
iam^ = max 



d t c PO (h 1 ,h 2 ) 



(12) 



with rp P ° < diam t ^ po < 2r t ( S po that make these mea- 
sures be equivalent. (Fig. [5]) plots in a color map the 
distance of each street to the topological center of Amiens 
that ends up to be a part of its highway-belt. This map 
gives a hierarchical vision of the space. There is no radial 
component of the increase of the topological distance: a 
scale of long streets serves the whole city, allowing the 
variation of the topological distance to be mainly local. 
Added to that the topological radius of the city grows 
very slowly with the size of the city (14 in the center of 
Amiens, 18 in the whole city that is eight times bigger). 
The topological efficiency l/dp P ° ' e " 9 (/io) is defined for 
each street and we can consider it is defined on each 
point of the city, beeing constant almost everywhere in a 
street. It defines a new centrality |llj on the network. 



C. Anisotropy 

In order to grant an efficient access to physical re- 
sources, the street system locally tends to be perpen- 
dicular to structuring elements as for example rivers or 
older streets. Then a city is not "isotropic". 
Let uq G R 2 be an arbitrary vector, taken as an angular 
origin. For angle a G [0, 7r] , 



**(a) = nc (c G C, Z(c,u ) € [0,a]) 



(13) 



where .) is the angle measure between two vectors in 
[0, tt]. It is a measure of the "total length" of streets in C 
that are oriented in directions [0, a]; in the special case of 
straight graphs, the impact of each street-segment is pro- 
portional to its length. From this, we define the angular 
density by \l/(a) = ^^"^ which representation describes 
the anisotropy of the city graph (Fig. [6j). We notice a 
fuzzy symmetry around the first bisectrix as a result of 
the streets local perpendicularity. For an isotropic city, 
the angular density would be a continuous and uni- 
form density ^i(a) = -. 

It would be useful to sum up this angular density as a 
single normalized indicator. We looked for a bound dis- 
tance measure between ^ and Wj. Since the observed 
distribution >F is discrete because of the limited num- 
ber of streets segments, measures like J \f — ~\ n highly 
depends on the size of the bean chosen to estimate the 
integral. 

Since the angle is defined modulo tt, we can "fold" 4" : 
*( 2 )(0) = ^{6)e i29 - / ^(u)e i2u du The inertia matrix 



h£H 



/<Re(*( 2 )pm(*( 2 >) 



J«e(*( 2 ))3m(*( 2 )) 
pm(*( 2 )) 2 



6 




FIG. 5. (Color online). The "topological maps" of Amiens (left) and its center (right). In each map, the red street is the 
topological center and the color of each street refers to its distance to that center. The maximum distance to the center is the 
radius of the map (18 in the whole city and 16 in the center). It is striking that this radius increases slowly with the size of 
the considered street system. We can even infer that the construction of surrounding highway belt (the found center) is made 
precisely to keep the topological radius of the city small. 




FIG. 6. Angular distribution of Amiens (top left : whole city - 
top right : center). We notice on both distributions a four- fold 
symmetry. The doubled angular distributions (bottom) look 
like ellipsoids which allows to define the anisotropy coefficient 
as the ratio of the eigenvalues of the inertia matrix: Ani = 
0.42 in the whole city and Ani — 0.71 in the center. 



is symetric and positive (from Cauchy-Schwartz's in- 
equality) with two eigenvalues Ai > A2 such that 



correctly defines an anisotropy coefficient. As for 
Amiens, its anisotropy coefficients varies from 0.42 in the 
whole city to 0.71 in the center. 



D. Street length 

The second order topology leads to think that a city 
organizes into a hierarchical way with long and no so long 
streets. Consequently we do not expect an exponential 
decay in the distribution of street lengths. In Amiens, the 
distribution of street lengthes L (Fig. [7]) is well-fitted by 
a mixture of Log-Normal laws: 

logi ~ p_.7V(m_,cr_) + (1 -p_)N(m + ,<y+) (15) 

with m_ < m+. The identification of this model is per- 
formed with an Expectation Maximization algorithm. 
The log-scaling reveals that there is no evident length 
scale in a city, there is no pre-existent typical street. A 
city auto-scales, its dynamics is purely multiplicative and 
could be the result of new streets cutting through former 
blocks or extending space at the exterior of the city. 
As for the bimodality, a town-planning explanation is 
that several transportation modes follow each other 
through time and their superposition creates modes in 
the distribution of L. 



Ani = 1 — 



(14) 
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Angoulme 


0.80 


0.28 


0.46 


300 


1628 


14 


0.21 


0.73 


0.37 


0.29 


39 


261 


7 


0.2 


Avignon 


0.85 


0.23 


0.59 


625 


3348 


13 


0.12 


0.82 


0.30 


0.73 


240 


1607 


9 


0.12 


Caen 


0.79 


0.29 


0.53 


485 


3045 


11 


0.08 


0.76 


0.34 


0.67 


128 


797 


9 


0.22 


Carcassonne 


0.86 


0.20 


0.66 


483 


1997 


20 


0.17 


0.72 


0.38 


0.89 


66 


296 


6 


0.42 


Dijon 


0.75 


0.33 


0.39 


558 


2605 


14 


0.14 


0.66 


0.42 


0.74 


149 


860 


8 


0.33 


Grenoble 


0.74 


0.32 


0.61 


361 


1638 


10 


0.17 


0.70 


0.40 


0.74 


119 


576 


6 


0.29 


Lyon 


0.66 


0.47 


0.53 


837 


3606 


15 


0.17 


0.53 


0.51 


0.93 


183 


606 


6 


0.24 


Rennes 


0.82 


0.26 


0.79 


625 


3538 


15 


0.20 


0.77 


0.30 


0.80 


192 


1041 


7 


0.17 


Rouen 


0.71 


0.38 


0.69 


348 


1770 


17 


0.19 


0.66 


0.43 


0.85 


141 


788 


7 


0.13 


Troyes 


0.81 


0.28 


0.87 


230 


1079 


9 


0.17 


0.67 


0.42 


0.92 


48 


248 


6 


0.41 



TABLE I. Features of 10 French cities and their centers. r n is the organic coefficient, M4 the meshedness coefficient, Ani the 
anisotropy, uc(C) the total street length, %H the number of streets in the hypergraph generated wit R w /io, r topo the topological 
radius of that hypergraph and res the root mean square between the street length distribution and its bi-lognormal fitting. 




Log of slroess len$h 



FIG. 7. The histogram represents the distribution of the log- 
arithm of streets length in the center of Amiens. The red plot 
is the fitting of this histogram with a gaussian mixture whose 
parameters are p~ = 0.5, m_ = 2.2, a_ = 0.3; m+ = 4.3, 
a+ = 1.2 . 

E. Synthesis: the city as a space division process 

(Tab. m summarizes the main indicators presented 
through this section for ten French cities and their ex- 
tracted centers. For French towns the organic coefficient 
rjv may seem surprisingly high, even for Lyon's center 
which is visually grid-like. The preponderance of degree 
3 intersections speaks in favor of seeing a city as the re- 
sult of a division process. The cutting of a block creates 
new intersections of degree 3 even if the block was formed 
of degree 4 vertices. 

The meshedness coefficient M4 is strongly correlated to 
rjv. In fact it is a second order refinement of rjy. We will 
show that it describes the shape of city graphs when rjv 
is constant equal to 1 in (Tab. 

The anisotropy coefficient Ani is the most discriminat- 
ing indicator: it ranges from 0.29 to 0.93. It also 



points out that the street system within a city is not 
homogenous: for instance j4mso(Lyon) = 0.53 and 
Aniso{Ljan. Center) = 0.93. 

The total street length /ic(C) and the number of streets 
jjif measure the size of a city. In fact the ratio between 
Hc{C) and provides an indication on the straightness 
of the city which is geometrical information. 
The topological radius r topo is very small compared to a 
measure of the city size: r topo <C $H. r topo may scale as 
the logarithm of §H, like a small-world complex network. 
This scaling may also be the consequence of a division 
process and/or of an upper system of long streets that 
covers the whole city's stretch constraining the topologi- 
cal distance's variations to be mainly local. However the 
topological radius is an indicator of the transportation 
performance of a given city: Carcassonne and Caen have 
almost the same total street length but the topological 
radius of Carcassonne is roughly twice Caen's. 
res the relative root mean square between the street 
length distributions and their bi-lognormal fitting is < 
0.2 with two exceptions (Carcassonne and Troyes's cen- 
ters) around 0.4. The lognormal scaling of street lengths 
is relevant in most of the cases. This also calls for a 
division (multiplicative) modelling of the city. But a di- 
vision process's length distribution would only have a 
single tail. 



In the following part (IV) we present a morphogenesis 
model for the organic city consisting in the duality be- 
tween division and extension of space. In (|v| we will 
show that this model reproduces the features we have 
presented above: a high organic coefficient, a small topo- 
logical radius and a local variation of the topological dis- 
tance plus a lognormal distribution for streets' lengths. 



IV. A STREETS BASED DYNAMICAL MODEL 

This section presents a model of the growth and de- 
velopment of a town. The town is reduced to its streets 
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and we build a dynamical model allowing adding street 
segments one after another. As in the previous parts, 
the spatial extension of the town and the geometry of 
the streets is of prim interest. As pointed out in [TJ], a 
city is above all an out of equilibrium system, that is to 
say a dynamical system observed at a random time of its 
development. 

Our model is first on three assumptions, two principles 
(installation and connection) and a few parameters. The 
whole giving a coherent and consistent vision of the prob- 
lem. We aim at building a model that can reproduce 
several limit cases of urban growth but also point out 
continuity between them. The principles and parameters 
we use are meaningful, expressed in an interface language 
that allows the mathematical and physical communities 
to exchange with town planners, architects and social sci- 
entists. 

This part develops the model in the quite general case 
of an organic development of the city on flat lands for 
which we can easily translate our assumptions into ana- 
lytical procedures. 



A. Hypothesis 

As a dynamical system and a geometrical graph, we 
will see a city as a function C : R+ — > Qc C Cd with 
C{t) — {V(t),E(t}}. Then we make the following postu- 
lates on the evolution of C : 

PI: A city is the result of a sequence of operations oc- 
curring at increasing times (ti)ieiN such that : 

C(t)=C(ti) Vte[t h t i+1 [ 

P2: Infrastructures are conserved: 

C(h)CC(t 2 ) if *i <f2 

P3: There exist two functions Pt (price) and V* (poten- 
tiality) such that the city is a compromise between 
them: 



C(t 



argmm 
c C C(t) 
P t (c-C(t)) <0 



V t (C(t),c) 



Functions C and P are not obvious to define. They 
should be in a "microscopic" point of view aggregation of 
economical parameters. We can avoid developing them 
if we observe a city's growth is determined by "macro- 
scopic" insights: 

Its planning: A city may be organic (the sum of local 
and independent phenomena: streets are added in- 
dependently with no visibility on a global planning) 
or centralized (a global authority decides of the co- 
herent and simultaneous addition of several streets 
on a large surface). 



Its construction: The capacity to add new elements to 
the map and to build new streets. 

Its organization: From a random settlement to a 
highly structured one. 

Its sprawling: A city has to make a compromise be- 
tween its inner development and outer growth. 

We will consider here the case of organic growth. 



B. Organic growth of a city 

The algorithm below simulates a city's growth within 
the individual settlement hypothesis. 
Under this assumption, each settlement (generic term to 
designate a commercial infrastructure, a private individ- 
ual...) is added at a given time and at a given location 
then connects to the existing infrastructures. 
The main idea is that the city C (t) induces a spatial po- 
tential field describing the attractiveness of any point of 
the available space. A new settlement (either an indi- 
vidual settler or a facility) has its own policy ( |IV B 3 ) of 
choice with respect to this potential. After having cho- 
sen its location, it connects to the existing street system. 
This model explicitly decouples the problems of position- 
ing and of connecting. 

In this section we will model how the geometry of the 
current city induces a potential field, how a new settle- 
ment can be connected to the city and at last how we can 
tune the behavior of settlements by a few parameters. 



1. Potential field 

For each point x in the available space the potential 
P(7_j. quantifies to what extent s is a good choice to 
locate a new center. This potential should mimic the 
following ideas: 

• A large scale behavior such that the global attrac- 
tion of a part of the city should be proportional 
to the global mass of infrastructures in place and 
slowly decrease with some distance d : Pc^x oc 

d~r{x,C) 

• A very short scale behavior that forbids a new 
center to be located on existing infrastructures: 

PC-yx = +00 

• A medium scale deduced from the two previous 
ones, that should display some local minimum. 

Thus among several possible fields whe choose (Fig. [8]): 



Pi 



C^rx — 



P 



due 



,(a?,C) ^[d 

'III I'll 

(x, C) 



\/dj_(c,x) 
(16) 
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FIG. 8. (Color online). The level lines of the potential field 
for a city reduced to a single segment of length lwith Ao = 1 
and /3 = 10. 



d m in(x,C) = min ce c' (x, c) is used in (Eq. 16) so that 
the rejecting zone is hard : there is a tube around the 
city where new settlements are impossible. The radius of 
this tube is Ao = (a/f3) 2 . 

dx(x,c) is the ||.||i norm in the local basis formed by 
the unitary tangent and the normal to C in c. The use 
of such a distance simplifies integral calculus compared 



to Euclidian distance. Towards oo, Pc^x 



d{x,C) 



and 



between those two extreme positions, interferences be- 
tween streets segments produce local minima. To choose 
parameters a and /3, one sets Ao the hard rejection ra- 
dius and /3 the long-range influence. The choice of /3 
influences the local geometry of the city but won't be 
discussed here. 

Among all possible potentials, we picked one that fulfils 
the conditions we set and that allows an explicit calculus 
of the integral. Further discussions would address the 
choice of the used distance and the decay exponent 7. 



close from each other, we introduce the relative neigh- 
borhood. In a general way, if P is a point and E a points 
set. Then s £ Els said to be in the relative neighborhood 
of P (s £ RN[E\\P}) if and only if Q3] 

Vu £ S, d{P, s) < max{d(w, P), d(u, s)} (19) 

ie there is no point both closer to s and to P. All candi- 
dats to become new connections are segments from x to 
RN{V x \ c \\x): [x,RN(C\\x)}. 



Parameters 



The tuning of those parameters will be discussed in 


a. Organization The global field induces minima. 
These minima represent points where it is the most in- 
teresting to settle. 

The question is to find a parameter P e that describes 
whether the city is organized or not. The idea is that 
when a city is organized it sticks strictly to optimal set- 
tlement places and when it is purely unorganized, new 
settlements are added at random without any influence 
of the potential field. 

Then a new settlement is selected by a Monte Carlo 
method with a number n of iterations and the new point 
is be chosen as X = Argmin n P(Xi). For random cities 
n is close to 1 and for organized one it is much higher. 
Let W be the area of a part of the plan that con- 
tains the current city, let Xi,...,X n be n points on 
this part, uniformly and independently chosen. And let 
X = argmin P(JQ). Then let 



P e = P(\X - Argmin P)\ < e) 



(20) 



represents the probability that the Monte Carlo method 
throws a point in a radius e of a local minimum. 
We want to give P e as an input parameter and traduce it 
into an iterations number. Let N be the number of local 
minima. If e is quite small: 



2. Connection 



Pr 



W - N.ir.e' 
W 



(21) 



Once a settlement is added in a location x, it links to 
the existing network C. Not all connections are eligible. 
From a point x we define the visible set of points: 

V xlc {x&C,[cx]nC = {c}} (17) 

And the optimal set of points from x of a part E of C: 

E x = {e £ E 3e | We' £ E n c ® B e , d(x, e) < d(x, e')} 

. (18) 

New connections are made between x and points in the 
optimal visible set V x \c- This one being a finite set in- 
cluded in C U (xJ-C) where X-LC is the set of orthogonal 
projections of x on the city. To avoid connections too 



U - 2 ( 22 ) 

e z iv.7r 

n is estimated roughly by noticing that a local minimum 
is often due to the interaction between two close streets 
segments: n ~ Z.v since 3 is roughly the average connec- 
tivity number of intersections. 

b. Connection and construction There are typically 
about four or five streets segments in [x, RN(V+\\x)] for 
a new center x. If the city shapes as a slum it would 
be tree like so we link the number of streets segments 
indeed added with the construction to £ [0, 1] of the city. 
We sort segments in [x, RN(V + \\x)] by increasing length : 
(sx, s n ). si is drawn with probability 1. n' ~ B(uj,n — 
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1) + 1 and segments S2,~s n < are also added. If uj = 1 
every admissible segment is added and if uj = only the 
shortest one. 

c. Sprawling When constructing with a rejection ra- 
dius Ao, the city gets a typical mesh width. If at a par- 
ticular urban operation a potential field with a rejection 
radius of KAq with K > 1 is considered then the city's 
inner meshes will appear as filled up with the rejection 
zone of this potential field and new points of interest will 
position outside of the city. 

With this observation we will consider that in a propor- 
tion f ext centers are added with respect to a potential of 
rejecting radius K ext X . This creates foils at the outskirt 
of the city and thus an extension of the city that repre- 
sent for instance an industrial zone which need a large 
surface. 

d. Refinements To enhance the realism of this 
model, some empirical parameters are added. 

The length of a new street segment is bounded to l max = 
fc; mai Ao- This can avoid too long costly connections (pos- 
sibly Z max = oo). 

The windows W out of where new settlements are chosen 
can be whether definitely set (Fig. |9| whether dynami- 
cally change with the overall city (Fig. 11 ) and (Fig. 13 1. 
Since a Monte-Carlo method is used to pick new centers, 
there is a very little probability that a new settlement 
is added in line with an existing street. If geometrically 
this does not have many consequences, it may strongly 
change the local topology. That is why, if an orthogonal 
projection is in a radius c\q with typically c ~ 0.3 of 
an intersection that is visible from the center then this 
orthogonal projection is removed. This rises the vertices 
degree and allows longer streets. 



V. A FEW SIMULATIONS 

To summarize an individual simulation of the city 
growth we need to provide our algorithm with several 
parameters: 

1. The number of settlements: N 

2. The organization probability P e 

3. The radius of the rejecting tube: Ao 

4. The long scale influence: ft 

5. The construction: to 

6. The sprawling factor K ext and the sprawling prob- 
ability f ext 

Notice that only four parameters will actually shape the 
simulated city (P e , ft, uj, f ext ) the others being scaling 
parameters and that the influence of ft won't be discussed 
here. 



A. Simulations with constant parameters 




FIG. 9. Simulations of the morphogenesis model with con- 
stant parameters with variation of the organization P e and 
the construction uj. On each thumbnail, the rejecting radius 
is Ao = 10m, there is no sprawling: k = 0, the number of 
settlements is N — 80 and the available space is bound in a 
square with sides of 400m. The red and bold segment repre- 
sents an initial street-segment and a scale of 20 meters. 



(Fig. [9]) shows the result of 16 simulations. The organi- 
zation probability P e and the construction uj are varying 
jointly when the same number of operations N = 80, the 
same rejecting radius Ao = 10m, the same available space 
( a square with an area of 1.6 km 2 ), the same initial city 
(a segment of length 20 m at the center of the available 
space) and the same extension probability of O are used. 
The first result is that this model is able to reproduce 
very different type of growth with very few "physical" 
parameters. 

We observe on this matrix representation that the 
meshedness M4 (|III| is an increasing function of both 
P e and uj (Tab. |II[). This result has been obtained by av- 
eraging the meshedness coefficient of 30 simulations for 
each couple (P e ,uj). Added to that the standard devi- 
ation of M4 for each couple is of 4 percent so that this 
coefficient is characteristic of the conditions of simula- 
tion. Contrary to that, the anisotropy coefficient Ani is 
almost the same in each case (between 0.31 and 0.46) 
with a large standard deviation of 20 percent. This Ani 
is quite large in the absolute: for the first iterations some 
directions have to be arbitrary chosen, which creates fa- 
vored directions. But it is the same order of weight as for 
the most isotropic French towns. Of course the organic 
ratio rpf is in every case close to 1. 

When uj ~ 0, the resulting simulations are to be com- 
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TABLE II. Variation of the meshedness coefficient M4 = 2(e— 
v + l)/(2.u.(l - N(2)) - 5) for the 16 simulated cities. Each 
result is the averaging of 30 simulations. The variation for 
each case is almost constant equal to 0.04. M4 is a increasing 
function of both P e and ui. 



pared to the Saffman- Taylor instability. It seems when 
uj ~ that only a bounded number of ramifications are 
possible from the initial segment (4 on this figure) as 
if first created branches shielded the initial center from 
newer ones. When u > 0, the resulting cities are to 
be compared to crack patterns: their dynamics follows a 
logic of division / subdivision of space. Interestingly, the 
Chinese town of Xi'an (Fig. 



10 ) has grown on a regular 



111 presents 
K f = 0.1, 



grid with a large mesh length with several populations 
that have different characteristics. The result is a gradi- 
ent of meshedness coefficient, from almost in the south- 
west to almost 1 in the top- right. (Fig. 
a city evolving with A = 10m, K = 10, 
P e = 0.8, oj = 0.7. The regular need of larger surface 
for activities such as industries, big institutions, etc. is 
well reproduced here. During the history, as the develop- 
ment of the city center progresses, it eventually absorbs 
the peripheral larger surfaces, to split them into smaller 
surfaces, with thus new larger places appearing at the 
new periphery. This reproduces and explains the situa- 
tion of economical zones always outside at the periphery 
of towns. It explains as well the successive subdivision of 
space, that leaves so many traces, first in the log normal 
distribution of streets length but also in the hierarchical 
distributions of streets (Fig. 12 1. For this simulation, the 
ratio r/v is equal to 0.93 so that the term "organic" fits. 
The meshedness coefficient M4 = 0.48 is quite close to 
Amien's (between 0.41 and 0.54) as the anisotropy (0.69 
to be compared to 0.71 in the center of Amiens). 
(Fig. 12 1 shows that the morphogenesis reproduces the 



small topological radius and the log-scaling of street 



lengths observed in real cities (HIE). We have run 20 
simulations with the same parameters. The root mean 
square distance between the length distribution and its 
best bi log-normal fitting range between 0.19 and 0.26. 
This compared to the result of the fittings for real French 
towns (less than 0.2 when the model is good and more 
than 0.4 when it is false) permits to claim that the model 
reproduces the street length scaling. 




FIG. 10. The Chinese town of Xi'an in 1949, whose various 
subdivision patterns inside a regular grid recall variations in 
the parameters of the morphogenesis model. 




FIG. 11. (a, b, c, d) : four steps in the development of 
the city (e )with 600 urban operations. For this simulation, 
A = 10m, K = 10, K f = 0.1, P e = 0.8, lo = 0.7 and the 
windows is adapted to the size of the current city. . The 
main phenomenon at work is the dynamics between the inner 
development and the extension of the city that creates two 
hierarchical scales. 



B. Simulations with varying parameters: the city's 
history 

Constant parameters are not realistic to model a real 
city, this one being shaped by its history which is from 
the morphogenetic point of view a variation of input pa- 
rameters. 

We represent the history of a city by a piecewise constant 
function t (P e ,w,? max , f extl K, fj, A ). 
For instance, the city of (Fig. |l~3"| ) has been obtained by 
simulating at first a city with a low construction and no 
sprawling [uj — 0.2 and f ext = 0) and then changed to a 
sprawling and constructed city (a; = 0.8 and f ex t — 0.15). 
The simulation starts with two perpendicular streets with 
a length 20 times larger than A = 10m. These pre- 
existing streets are structuring elements as could be a 
river. This kind of variation in parameters recalls Kas- 
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FIG. 12. (Color online). The distribution of the logarithm of 
streets length (top) and the second order topology (bottom) 
for the resulting city of simulation (Fig. 111. As for real 
cities (Amiens), the second order topology presents a bounded 
hierarchical representation of the city and the streets length 
is well- fitted by a mixture of log-normal random variables. 



bah in Morocco where the historical center of the city is a 
souk. This model is the first step of a very simple model. 
We can see that with only one type of event (new settle- 
ment) and a few parameters a great variety of structures 
can be obtained. More refinements can be added. The 
first one is the planned creation of a highway belt as in 
Amiens or an enclosure using the punctual addition of 
the convex hull of the current city. The second possibil- 
ity would be to add also the "Hausmann" effect, allow- 
ing splitting pre-existing streets with new street patterns. 
We could also consider distinct populations with interac- 
tion rules that build a city in the same time. These three 
points are going to be developed in a second version of 
the morphogenesis model. 




FIG. 13. A city developing with varying parameters from two 
perpendicular long structuring elements. The historic center 
of the city has been built with a low construction parameter 
and no sprawling (w = 0.2 and f ex t = 0) to recall the tree 
aspect of a central souk in the Kasbah. Then parameters 
are changed to ui — 0.8 and f ex t = 0.15 which products an 
industrial crown. 



VI. CONCLUSION 

We have reduced cities to the map of their street seg- 
ments and shown that a lot of information can be deduced 
from this representation without additional data such as 
population dispersion, width of streets, ground specific 
use. For instance we can find back streets, characterize 
topology and shape or define a centrality. To this we 
have introduced the notion of geometrical and straight 
graph with a canonical hypergraph structure to define 
difference between streets segments and streets. A mea- 
sure allows seeing a city's map as a "continuous graph" 
or likewise as an object both relational and geometrical 
0. 

While [3 considered a city as a pure graph and [IT] took 
into account its spatiality, we have explored the geometri- 
cal aspect of the city, its topology being only the skeleton 
that holds it up. 

From this point of view we have shown that despite an 
evident diversity on their overall shapes (anisotropy first- 
order topology), a few fundamental rules can explain 
cities' general morphogenesis (IV). The urban infrastruc- 



ture differentiates to adapt to the local geography and to 
fulfill the constraints of lodging people while maintain- 
ing their efficient transportation. The model developed 
in (IV) and illustrated in (M) is based on the division / 



extension of space principle observed for French towns in 
(III). It shows that structural properties of cities stand 



out of the local constraints and behaviors that define the 
dynamics of cities. For instance the log scaling of result- 
ing street system is a global, non trivial property that 
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validates the model as well as the small observed topo- 
logical radius. Even in the organic case when there is 
no global and coherent plan the topological radius of the 
city increases slower than the size of the city. 
To simulate cities' dynamics, we have uncoupled the 
space potential induced by the current infrastructure, the 
policy of connections and the freedom a new settlement 
(generic term to refer to a commercial infrastructure, a 
private individual...) ought to take on the two previous 
rules. 

Further than the interest of a physical modelling of their 
object of study, our work can be applied by town plan- 
ners and social scientists: the morphogenesis can help 
for semi-supervised planning and the analysis of cities al- 
lows detecting abnormalities on a map (the second-order 
topology representation is particularly effective). It can 
also be applied in engineering to test a technology that 
strongly depends on the urban infrastructures. 
This morphogenetic model calls for a lot of outcomes: 



study of several potential fields, comparison to a large 
data basis of existing cities, following of the evolution of 
parameters as the city grows. The important element is 
that morphogenesis implements a space division and ex- 
tension process. Since we used arbitrary tunnings and 
yet obtained realistic results, we have shown that this 
class of process is robust to model cities, whether it be 
in the centralized case or in the organic one. 
In the centralized case, the reason wants the street sys- 
tem to be the homotopy of a square grid to adapt to the 
local geography and infrastructure, which can be seen 
afterwards as the result of a correlated division process. 
Indeed in the organic case, the city's layout comes from 
the duality between the city's expansion and the grain- 
ing of former large cadastres by new settlements. And 
from this local division process emerge some non trivial 
global phenomena (log scaling of streets, low topological 
radius). One can thus say that the division of space is a 
natural response of cities to fulfill their functional goals. 
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